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ABSTRACT 

High-resolution N-body simulations are used to investigate systematic trends in the mass profiles and 
total masses of clusters as derived from 3 simple estimators: (1) the weak gravitational lensing shear 
field under the assumption of an isothermal cluster potential, (2) the dynamical mass obtained from 
the measured velocity dispersion under the assumption of an isothermal cluster potential, and (3) the 
classical virial estimator. The clusters used for the analysis were obtained from simulations of a standard 
cold dark matter universe at z = 0.5 and consist of order 2 x 10 5 to 3 x 10 5 particles of mass m p ~ 10 10 M Q . 
The clusters are not smooth and spherically symmetric but, rather, the mass distribution is triaxial and 
significant substructure exists within the clusters. Not surprisingly, the level of agreement between the 
cluster mass profiles obtained from the various estimators and the actual mass profiles is found to be 
scale-dependent . 

We define the total cluster mass to be the mass contained within a 3-dimcnsional radius T2oo of the 
cluster centers, where r2oo is the radius inside which the mean interior overdensity is equal to 200. Under 
this definition the classical virial estimator yields a good measurement of the total cluster mass, though 
it is systematically underestimated by ~ 10%. This result suggests that, at least in the limit of ideal 
data, the virial estimator is quite robust to deviations from pure spherical symmetry and the presence 
of substructure within a cluster. The dynamical mass estimate based upon a measurement of the cluster 
velocity dispersion and an assumption of an isothermal potential yields a poor measurement of the total 
cluster mass, underestimating it by ~ 40% for the case that a v is computed from an average over the 
entire cluster. The weak lensing estimate yields a very good measurement of the total cluster mass, 
provided the mean shear used to determine the equivalent cluster velocity dispersion is computed from 
an average of the lensing signal over the entire cluster (i.e. the mean shear is computed interior to a 
projected cluster radius of i?20o)- 

Subject headings: cosmology: theory — dark matter - gravitational lensing — large-scale structure of 
the Universe — galaxies: clusters: general — methods: numerical 



1. INTRODUCTION 

Rich clusters of galaxies constitute the largest 
gravitationally-bound objects in the universe and the his- 
tory of their formation is a potentially powerful test of 
the viability of differing large-scale structure models. In 
particular, the time evolution of the cluster mass function 
is expected to be a good discriminator between low- and 
high-density universes, and, additionally, it provides a con- 
straint on the degree of bias between the galaxy and mass 
distribution (e.g. Bahcall et al., 1997; Fan et al., 1997). 
Owing to the lack of complete, uniform, mass-selected cat- 
alogs of galaxy clusters out to large depths (z ~ 0.5 to 
1.0), the evolution of the cluster mass function is not con- 
strained especially well at present. It is expected, however, 
that significant effort will soon be devoted to developing 
such catalogs, particularly using deep wide-field imaging 
of weak gravitational lensing of distant field galaxies by 
intervening clusters. The ability of these future catalogs 
to constrain models of structure formation rests heavily 
on the accuracy with which the cluster masses can be 
obtained and an understanding of any systematic biases 
present in the mass estimators themselves. 



The virial mass estimator is the method which has 
the longest history of application to galaxy clusters and 
it yields consistent results for cluster mass to light ra- 
tios in the range of 200h{M/L) Q to 4Q0h(M/L) Q (e.g. 
Zwicky 1933, 1937; Smith 1936; Schwarzschild 1954; Gott 
& Turner 1976; Gunn 1978; Ramella, Geller & Huchra 
1989; David, Jones & Forman 1995; Bahcall, Lubin & Dor- 
man 1995; Carlberg et al. 1996; Carlberg, Yee & Ellingson 
1997). While there are legitimate concerns that large clus- 
ters are not fully virialized, Carlberg et al. (1997a) have 
presented spectroscopic evidence which strongly suggests 
that the clusters in the CNOC survey are in equilibrium 
and, therefore, that the masses obtained within the virial 
radius are reliable. 

Nevertheless, because of its potential power to map the 
dark mass distribution within a cluster independent of the 
cluster's dynamical state, recently a great deal of effort 
has been devoted to measurements of the gravitational 
potentials of clusters via observations of weak lensing (e.g. 
Tyson, Wenk & Valdes 1990; Bonnet et al. 1994; Dahle, 
Maddox & Lilje 1994; Fahlman et al. 1994; Mellier et al. 
1994; Smail et al. 1994, 1995, 1997; Tyson & Fischer 1995; 
Smail & Dickinson 1995; Kneib et al 1996; Seitz et al. 



1 



2 



A Comparison of Simple Mass Estimators for Galaxy Clusters 



1996; Squires et al. 1996ab; Bower & Smail 1997; Fischer 
et al. 1997; Fischer & Tyson 1997; Luppino & Kaiser 1997). 
Relatively few clusters have been studied in detail, but a 
consistent picture of the dark mass distribution appears to 
be emerging from the weak lensing investigations. In par- 
ticular, the center of mass corresponds well with the center 
of the optical light distribution and the smoothed light dis- 
tribution traces the dark mass well. The lensing-derived 
mass to light ratios vary from cluster to cluster but bracket 
a broad range of 200h(M/L) Q to 8OO/i(M/L) , with most 
of the clusters falling in the middle of the range. 

The consistency of cluster masses obtained from inde- 
pendent methods such as lensing, virial analyses, or X- 
ray data (assuming pressure supported hydrostatic equi- 
librium) is very much in debate at this time. Cluster mass 
estimates obtained from observations of strong lensing can 
often exceed the X-ray mass by a factor of 2 to 3 (e.g. 
Miralda-Escude & Babul 1995). This particular discrep- 
ancy is likely due to the failure of the assumption of hydro- 
static equilibrium and Waxman & Miralda-Escude (1995) 
have thus proposed the existence of multiphase cooling 
flows in the centers of rich clusters. Based on analyses of 
weak lensing over considerably larger cluster radii, how- 
ever, some studies conclude that there is quite good agree- 
ment between the lensing and X-ray masses (e.g. Squires et 
al. 1996ab; Smail et al. 1997) while others claim a signifi- 
cant disagreement in which the lensing mass systematically 
exceeds the X-ray mass (e.g. Fischer et al. 1997; Fischer & 
Tyson 1997). Additionally, the weak lensing mass estimate 
is often found to exceed the mass obtained from the virial 
estimator by a factor of order 2 (e.g. Fahlman et al. 1994; 
Carlberg et al. 1994; Smail et al. 1997) but in some cases 
good agreement between these two independently derived 
masses is found (e.g. Fischer et al. 1997). 

One of the troubles associated with a fair assessment of 
independent cluster mass estimates is that different tech- 
niques tend, out of necessity, to be applied at different 
cluster radii (i.e. it is not always possible to investigate the 
gravitational potential over the entire cluster via one par- 
ticular estimator simply due to lack of data on appropri- 
ate scales). Systematic biases inherent in any given mass 
estimator may well be scale-dependent and amongst dif- 
ferent estimators the form of the scale dependence is likely 
to vary. Therefore, a considerable effort will be required 
in order to reconcile all of the outstanding discrepancies 
amongst independent cluster mass estimates. 

Few observational investigations have been able to place 
direct constraints on the radial mass profiles of clusters 
to date. Bonnet et al. (1994) detected a coherent weak 
lensing shear due to the cluster C10024+1654 out to a ra- 
dius of r ~ 1.5/i Mpc and from their observations they 
showed that the underlying cluster mass profile was con- 
sistent both with an isothermal profile and a steeper de 
Vaucouleurs profile. Similarly, Fischer & Tyson (1997) 
found that the weak lensing shear field of RXJ 1347.5- 
1145 yielded a density profile that was consistent with 
isothermal. Tyson & Fisher (1995), however, found that 
the density profile implied by weak lensing observations 
of A1689 was steeper than isothermal on large scales 
(200/i _1 kpc £ r £ lft^Mpc). Similarly, Squires et al. 
(1996b) found that the density profile of A2390 implied 
by the weak lensing shear field was consistent with isother- 
mal on small scales (r < 250/i -1 kpc), but on larger scales 



was better described by a profile steeper than isothermal. 
In contrast to the weak lensing results, however, Carlberg 
et al. (1997b) found that the velocity dispersion profiles of 
the CNOC clusters gave rise to a mean cluster mass profile 
that was fit very well by a Navarro, Frcnk & White profile, 
i.e. shallower than isothermal at small radii and isother- 
mal at large radii (e.g. Navarro, Frenk & White 1995, 1996, 
1997). 

In this investigation we use high-resolution N-body sim- 
ulations of rich clusters to investigate systematic trends 
in both the total cluster masses and the radial mass pro- 
files as derived from three simple estimators: (1) the weak 
lensing shear field under the assumption of an isothermal 
cluster potential, (2) the dynamical mass obtained from 
the measured velocity dispersion under the assumption of 
an isothermal cluster potential, and (3) the classical virial 
estimator. The simulated clusters are very massive and 
thus do not constitute an average, unbiased sample of ob- 
jects. They do, however, correspond to the largest clus- 
ters likely to form in a standard CDM universe and are 
objects which would certainly be detectable as weak grav- 
itational lenses. The N-body simulations of the clusters 
are discussed in §2 and the weak lensing and dynamical 
properties of the clusters are discussed in §3 together with 
the mass profiles obtained from the three estimators. A 
discussion of the results is presented in §4. 

2. THE NUMERICAL CLUSTERS 

The Hierarchical Particle-Mesh (HPM) N-body code 
written by J. V. Villumsen (Villumscn 1989) was used to 
simulate the formation of three rich clusters. The HPM 
code allows small- volume particle-mesh simulations to be 
nested sclf-consistently within large- volume particle-mesh 
simulations and by successively nesting many simulations 
within each other it is possible to obtain extremely high 
resolution in both mass and length within a small, local- 
ized region of a large computational volume (a "power 
zoom" effect). The code is, therefore, especially useful for 
the simulation of the formation of objects such as individ- 
ual clusters. In particular, using the HPM code to simu- 
late the formation of clusters obviates the need for compu- 
tations which utilize "constrained initial conditions" and 
those which simulate at high resolution the evolution of 
density peaks that have been excised from the initial condi- 
tions of a large computational volume. While the largest- 
volume, lowest-resolution grid in an HPM simulation uses 
periodic boundary conditions, the smaller- volume, higher- 
resolution grids use isolated boundary conditions, allow- 
ing mass to flow in and out of the higher-resolution grids 
over the course of the simulation. Therefore, unlike con- 
strained initial conditions or peak-excision simulations, 
clusters simulated with HPM are guaranteed to accrete 
all of the mass that they should accrete if one simply ran 
a single large-volume simulation at a level of resolution 
comparable to that of the highest-resolution HPM grid. 

The HPM code uses a standard cloud-in-cell (CIC) inter- 
polation scheme, which results in an approximately Gaus- 
sian smoothing of the power spectrum with a smoothing 
length of r s — 0.8 grid cell (see, e.g., §6.6 of Blandford 
et al. 1991). The gravitational force is, therefore, softer 
than Newtonian on small scales but becomes Newtonian 
for length scales greater than or of order 2 grid cells (Vil- 
lumscn 1989). Due to the force softening we therefore 
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restrict our analyses to length scales greater than 2 grid 
cells. 

The clusters used for the present analysis are discussed 
in detail in Brainerd, Goldberg & Villumsen (1998) and 
here we present only a summary of the simulations in- 
volved. A standard Cold Dark Matter model (Qq = 1, 
A = 0, and H = 50 km/s/Mpc) was adopted and the 
present epoch (i.e. redshift, z, of 0) was taken to corre- 
spond to er 8 = 1 where 



cr s = 



P 



(1) 



This is a model which is somewhat under-normalized com- 
pared to the COBE observations (e.g. Bunn & White 1997) 
and over-normalized compared to the abundance of rich 
clusters (e.g. Bahcall & Cen 1993; White, Efstathiou & 
Frenk 1993; Eke, Cole & Frenk 1996; Viana & Liddle 
1996). The simulations began at er 8 = 0.033 (correspond- 
ing to a redshift of 29) and were evolved forward in time 
to cr 8 = 1.0. 

The formation of the three most massive (i.e. "richest") 
clusters contained within a single cubical volume of comov- 
ing side length L = 400 Mpc was followed at high resolu- 
tion. The large, primary computational volume common 
to all three clusters was a standard particle-mesh simula- 
tion consisting of 256 3 grid cells and 128 3 particles. Each 
cluster in turn was simulated at increased resolution by 
nesting two smaller, higher-resolution grids successively 
within the large, primary simulation volume. These two 
small, higher-resolution grids were centered on the cen- 
ter of mass of the particular cluster being simulated, used 
256 3 grid cells each and had comoving side lengths of 
L = 66.6 Mpc and L = 16.7 Mpc, respectively. Since 
isolated boundary conditions are used in the small- volume 
simulations, the number of particles physically inside the 
small grids varied over the course of the simulations. 

The particle mass in the large computational volume 
common to all three clusters was m p = 2.1 x 10 12 M Q , 
while for the smaller grids unique to the individual clus- 
ters, the particle masses were m p — 7.8 x 1O 1O M and 
m p = 9.8 x 10 9 M Q in the grids with L = 66.6 Mpc and 
L = 16.7 Mpc, respectively. Dynamic ranges of <~ 4.5 x 10 8 
in mass and ~ 6000 in length were thus achieved by nesting 
the smaller simulations within the primary computational 
volume. 

Throughout the present analysis we shall define the clus- 
ters to consist of all particles within a radius r 2 oo of the 
centers of mass, where r 20 o is the radius inside which the 
mean interior mass overdensity is 200: 



5p 



(r 2 oo) = 200 



(2) 



(see, e.g., Navarro, Frenk & White 1997, 1996, 1995). The 
cluster mass estimates were computed using only parti- 
cles from the three highest-resolution grids (i.e. the grids 
with L = 16.7 Mpc and m p = 9.8 x 1O 9 M ) at an epoch 
corresponding to a redshift of 0.5 (os = 0.67 for our nor- 
malization) . All of the particles within r 2 oo of the cluster 
centers were excised from the highest-resolution grid and 
were then used to compute the mass profiles. The to- 
tal number of particles per cluster located within r 20 o at 



z = 0.5 are: 192346 ("cluster 1"; r 200 = 2.1 Mpc, proper 
radius), 288641 ("cluster 2"; r 200 = 2.4 Mpc, proper ra- 
dius), and 310310 ("cluster 3"; r 20 o = 2.5 Mpc, proper 
radius). Within r 20 o the clusters contain a significant 
amount of substructure and have median projected ellip- 
ticities of 0.3. Cluster 1 is nearly prolate while clusters 2 
and 3 are nearly oblate. Although the density profiles of 
the clusters are fit well by Navarro, Frenk, & White (1997, 
1996, 1995) profiles, the values of the best-fit concentra- 
tion parameters obtained for the clusters are a factor of 
order 2 lower than the values predicted by the Navarro, 
Frenk & White formalism for objects in the identical mass 
range. For a full discussion of the above cluster properties 
see Brainerd, Goldberg & Villumsen (1998). 

3. RESULTS 

Two-dimensional projections of the clusters are shown 
in the top panels of Figs. 1, 2, and 3. The color scale 
shows the logarithm of the surface mass density (in units 
of M Q /kpc 2 ) and distances are given in proper coordi- 
nates for z = 0.5. The actual mass profiles of the clus- 
ters are shown in Fig. 4, where the top panel shows the 
mean 2-dimensional projected mass profile obtained from 
10 random projections of each cluster and the bottom 
panel shows the 3-dimensional mass profile of each clus- 
ter. Throughout we shall adhere to notation in which R 
refers to a proper radius projected on the sky and r refers 
to a 3-dimensional proper radius. In this notation, then, 
M(R) is the projected mass interior to a radius R on the 
sky and M(r) is the mass interior to a sphere of radius r. 
As expected from the work by Navarro, Frenk & White 
(1997, 1996, 1995) on the relatively generic shapes of the 
density profiles of objects formed by dissipationless col- 
lapse, the mass profiles of the numerical clusters are not 
fit well by single power laws (see also Dubinski & Carlberg 
1991, Cole & Lacey 1996, and Tormen, Bouchet & White 
1997). Rather, a gently changing slope is observed, with 
the density profiles becoming roughly isothermal on large 
scales (> 1 Mpc in the case of our clusters). 



Figs. 1 , 2 & 3 are very large and have been omitted. 



Fig. 1: Top panel: the logarithm of the surface mass density of cluster 
1 as observed from a randomly-chosen line of sight. The units of the 
surface mass density are MQ/kpc 2 . The cluster consists of all particles 
in the highest-resolution subgrid that arc located within a radius r2oo 
of the center of mass. Bottom panel: the gravitational lensing shear, 7, 
obtained for the projected mass density shown in the top panel. The 
cluster was placed at a redshift of 0.5 and the shear that would be in- 
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duccd in a plane of sources at z — 1.0 was computed by tracing a regular 
grid of 4 X 10 light rays through the cluster. The color scale indicates 
the local value of log 10 7 while the orientation of the sticks indicates the 
orientation of the shear, <p. For clarity, the mean orientation of the local 
shear is shown on a coarse 10 X 10 grid. The angular scale of the figure is 
of order 11 X 11 . Figs. 2 and 3 are the same as Fig. 1, but for clusters 
2 and 3, respectively. 

In the following subsections we will compare the true 
mass profiles of the clusters to the mass profiles obtained 
from the three estimators. All of the mass estimators as- 
sume the clusters to be spherically symmetric and, addi- 
tionally, both the weak lensing and "isothermal" dynam- 
ical mass estimates assume that the cluster potential is 
approximately isothermal. 
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Fig. 4: The mass profiles of the clusters as computed directly from the 
distribution of particles in the highest-resolution grids. Top panel: the 
mean projected mass profile, computed from 10 random projections of 
each cluster. Bottom panel: the full 3-dimcnsional mass profile. The 
clusters arc roughly isothermal on scales ^ 1 Mpc. 

3.1. Weak Lensing Shear 

Observations of gravitational lensing provide potentially 
powerful constraints on both the total mass and the mass 
distribution within clusters of galaxies. The gravitational 
potential of the cluster systematically deforms the shapes 
of distant source galaxies that are seen through the lens- 
ing cluster. The result is a net ellipticity induced in the 
images of lensed galaxies and a net tangential alignment 
of the lensed images relative to the center of the cluster 



potential. 

Provided the distance traveled by the light ray is very 
much greater than the scale size of the lens itself, it is 
valid to adopt the "thin lens approximation" in order to 
describe a gravitational lens. Consider a lens with an arbi- 
trary 3-dimcnsional potential, $. In the thin lens approx- 
imation a conveniently scaled 2-dimcnsional potential, tp, 
is adopted for the lens (i.e. -0 is a scaled representation 
of the 3-dimcnsional potential of the lens integrated along 
the optic axis): 



Dds 
D d D s 



z) dz. 



(3) 



Here 8 is the location of the lensed image on the sky rela- 
tive to the optic axis, Dd s is the angular diameter distance 
between the lens (the "deflector") and the source, is 
the angular diameter distance between the observer and 
the lens, and D s is the angular diameter distance between 
the observer and the source. 

Having adopted this 2-dimensional lens potential, then, 
it is straightforward to relate the potential of the lens 
(through second derivatives of ip) directly to the two fun- 
damental quantities which characterize the lens: the con- 
vergence (k) and the shear (7). The convergence, which 
describes the isotropic focusing of light rays, is given by: 



. . d 2 ip d 2 ip 



(4) 



The shear describes the tidal gravitational forces acting 
across a bundle of light rays and, therefore, the shear has 
both a magnitude, 7 = \J^\ + and an orientation, (p. 
In terms of ip, the components of the shear are given by: 



7i (?) 



d 2 ^ 



del 



12(8) = 



d 2 ip 



39^62 de 2 



7((9)cos 2<p(0) (5) 



= 7(0) sin 2ip(0) (6) 



(e.g. Schneider, Ehlers & Falco 1992). 

A great deal of work has been done in recent years to de- 
velop methods by which a map of the surface mass density 
of a cluster can be reconstructed from observations of the 
distortions induced in the images of background galaxies in 
the limit of weak gravitational lensing, for which k << 1 
and I7I << 1 (e.g. Kaiser & Squires 1993; Bartelmann 
1995; Kaiser 1995; Kaiser et al. 1995; Schneider 1995; 
Schneider & Seitz 1995; Seitz & Schneider 1995; Bartel- 
mann et al. 1996; Seitz & Schneider 1996; Squires & Kaiser 
1996; Seitz et al. 1998). It is not the intent of this paper 
to explore these detailed methods of cluster mass recon- 
struction. Rather, we will focus on a very simple weak 
lensing analysis technique that is sometimes used to gauge 
the total mass of a cluster contained within a given radius 
without fully reconstructing the underlying density profile. 

The method invokes an assumption that the cluster po- 
tential may be represented adequately by an isothermal 
sphere. The actual density potentials of the numerical 
clusters are better represented by Navarro, Frenk & White 
profiles (see Braincrd, Goldberg & Villumsen 1998) than 
by singular isothermal spheres and, given the apparent 
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generality of the NFW profile, it is more likely that an 
NFW profile will better represent an actual galaxy clus- 
ter than will an isothermal sphere. Here we choose to 
adopt the isothermal sphere approximation for the analy- 
sis because this is the simplifying assumption that is most 
commonly invoked in the literature when cluster masses 
are estimated from observations of weak lensing without 
a full reconstruction of the density profile (see the refer- 
ences listed below). Here our goal is simply to quantify 
systematic effects due to the assumption of an underly- 
ing isothermal potential when the true potential is better 
approximated by that of an NFW-type object. 

An isothermal sphere is uniquely specified by a single 
quantity, the velocity dispersion (a v ), and the mass of an 
isothermal sphere contained within a 3-dimensional radius 
r is given by 

2<r 2 r 

M(r) = -g- (7) 

where G is Newton's constant. The total mass of an 
isothermal sphere within a radius R projected on the sky 
is given by 

M(R) = ^ (8) 

(e.g. Binney & Tremaine 1987). 

Since it is spherically symmetric, the isothermal sphere 
gives rise to a gravitational lensing shear field which is nec- 
essarily circularly symmetric and, in particular, the shear 
as a function of angular radius, 9, is given by 



7(0) 





'D ds ' 




[Ds\ 



(9) 



where c is the velocity of light and a v is the velocity dis- 
persion of the lens (e.g. Schneider, Ehlers & Falco 1992). If 
we consider an annulus of inner radius # m i n and outer ra- 
dius ^max centered on the center of mass of the isothermal 
sphere, the mean shear inside the annulus is given by 



7 



47T 



(ft, 



o- 1 . 



(10) 



That is, provided the cluster potential is sufficiently well- 
represented by an isothermal sphere it is possible to deduce 
its characteristic velocity dispersion directly from either a 
measurement of the shear at a given radius, r )(9) 1 or the 
mean value of the shear, 7, computed within some large 
large annulus. A measurement of a v by such a technique 
then leads to an estimate of the mass of the cluster within 
a given radius (e.g. Tyson, Wenk & Valdes 1990; Bonnet 
et al. 1994; Smail et al. 1994, 1997; Smail & Dickinson 
1995; Bower & Smail 1997; Fischer & Tyson 1997). 

In practice, an observed weak lensing shear only places 
a limit on the mass of cluster to within an additive con- 
stant (the so-called uniform density mass sheet degener- 
acy). The simple singular isothermal sphere mass estima- 
tor that we use here formally assumes that there is no 
such mass sheet present and that the observed weak lens- 
ing shear can be directly translated into a mass measure- 
ment via equations (7), (8), (9), and (10). Below we will, 
therefore, interpret the cluster shear fields in a manner 
consistent with the simple form of the mass estimator and 
we will not explicitly address the mass sheet degeneracy 



problem or its implications for an observed weak lensing 
shear. 

In this section we compute the shear fields of the nu- 
merical clusters and in §3.3 we will use these shear fields 
to investigate the systematic effects that the above weak 
lensing mass estimate has on the masses inferred for the 
numerical clusters. The shear fields of the clusters are 
determined directly by tracing regular Cartesian grids of 
2001 x 2001 light rays through the clusters. In the anal- 
ysis below we adopt the thin lens approximation and for 
a particular plane projection of a cluster we simply calcu- 
late the net deflection of each light ray due to all of the 
point masses contained within r 2 oo of the cluster center 
of mass. Note, however, that we ran a few test cases in 
which all particles inside a radius of 4 Mpc of the cluster 
centers were included in the ray trace analysis. The inclu- 
sion of the mass exterior to a radius of r 2 oo gave rise to a 
shear field interior to r 2 oo that was indistinguishable from 
the shear field obtained using only the particles interior 
to r2oo- That is, owing to the fact that the clusters are 
roughly axisymmetric and no large mass concentrations 
exist just outside the clusters, the shear interior to a pro- 
jected radius R is determined by the surface mass density 
interior to R. 

The clusters are located at a redshift of z = 0.5 and 
we consider a plane of sources at z = 1.0. (Although the 
redshift of the sources will affect the magnitude of the 
shear, it will not affect the velocity dispersion inferred in 
the isothermal sphere approximation and, therefore, the 
choice of source plane is essentially arbitrary for our anal- 
ysis.) The side lengths of the grids of light rays were taken 
to be L = 2r 2 oo so that throughout we compute only the 
shear interior to the virial radii of the clusters. At the 
redshift of the clusters, then, the side lengths of the grids 
correspond to and angular scale of order 11' x 11'. 

If we let the location of a light ray on the grid be given 
by (3 prior to lensing (i.e. (3 is the location of the light 

ray in the source plane) and we let 9 be the location of 
the light ray after having been lensed by all of the point 
masses (i.e. 9 is the location of the light ray in the image 
plane) , the components of the shear are then given by: 



(11) 
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(12) 



The 2001 x 2001 light rays define a grid of 2000 x 2000 cells 
and the shear at the centers of each of these cells can be 
determined from equations (11) and (12) above by finite 
differencing of the deflections of the four light rays which 
define the corners of the cell. 

The code used to compute the net deflections of the grid 
of light rays was tested by tracing the light rays through a 
number of singular isothermal spheres that were approx- 
imated by a set 250,000 point masses. The point masses 
were constrained to lie within a maximum projected radius 
of R = 2.7 Mpc and their masses were scaled appropriately 
so as to reproduce the correct values of M(R — 2.7 Mpc) 
for a set of isothermal spheres with values of a v in the 
range of 500 km/s to 1500 km/s. As with the simu- 
lated clusters, the isothermal sphere lenses were placed at 
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z = 0.5 and the source light rays emanated from z = 1.0. 
The net deflections of the light rays were evaluated and the 
radial dependence of the convergence, k(R), and shear, 
7(.R), was computed and compared to the analytic ex- 
pectations for infinite singular isothermal spheres having 
values of o~ v identical to the isothermal spheres that were 
approximated by the point masses. For the isothermal 
sphere we know k(R) = 'y(R), and in all cases good agree- 
ment was found between both k(R) and j(R) as computed 
individually from the ray tracing and between the ray trac- 
ing results and the analytic expectations (deviations < 1% 
of the analytic values). 

Shown in the bottom panels of Figs. 1, 2, and 3 are the 
shear fields corresponding to the 2-dimensional projections 
of the clusters shown in the top panels of these figures. 
The color scale shows the logarithm of the magnitude of 
the shear and the small sticks indicate its orientation. For 
clarity of the figure, we plot the mean orientation of the 
shear on a coarse 10 x 10 grid that was computed from 
an unweighted average of the local shear vectors obtained 
from the differencing of the displacements of the light rays 
(i.e. the sticks show a rebinning of the original 2000 x 2000 
grid of shear vectors onto a 10 x 10 grid). The visual agree- 
ment of the magnitude and orientation of the shear with 
the actual surface mass density of the clusters is as ex- 
pected; the shear is greatest in the densest regions of the 
clusters and is oriented roughly tangentially with respect 
to the cluster centers. The shear fields are not, however, 
circularly symmetric and reflect both the overall ellipticity 
of the clusters and the substructure within them. 
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Fig. 5: The mean gravitational lensing shear for the clusters as a func- 
tion of projected radius. Two-dimensional shear fields were determined 
for 10 random projections of each cluster, from which the average ra- 
dial value of the shear was computed in independent bins of radius R 



centered on the cluster center of mass. The error bars show the formal 
standard deviation in the mean between the 10 projections. For compar- 
ison the dotted line indicates the shape of the shear profile expected for 
an isothermal sphere lens, ■y(R) oc R~ 1 . 

Each cluster was viewed at 10 random orientations and 
a mean radial shear profile was computed from the full 
2000 x 2000 grid of shear vectors. The results are shown 
in Fig. 5, where the error bars indicate the formal standard 
deviation in the mean between the 10 random projections. 
Also shown for comparison is the radial shear profile ex- 
pected for an isothermal sphere (i.e. " (1?) oc R , cf. equa- 
tion 9 above). Below a scale of ~ 1 Mpc the radial shear 
profiles of the clusters behave as j(R) oc R~ 5 , while on 
larger scales the variation of 7 with R is roughly isother- 
mal, j(R) oc R . Given the mass profiles shown in Fig. 4, 
this is precisely the behavior we would anticipate for the 
shear profiles. This behavior will, however, cause system- 
atic errors in the cluster masses inferred from the shear 
fields under the assumption of isothermal cluster poten- 
tials. 

3.2. Velocity Dispersions 

Under the assumption of an isothermal cluster potential, 
the masses of the clusters can be determined from mea- 
surements of their velocity dispersions alone (e.g. equa- 
tions 7 and 8 above). The isothermal sphere is charac- 
terized by a single, constant value for the velocity dis- 
persion and in this section we investigate the degree to 
which the measured cluster velocity dispersions vary with 
distance from the cluster centers of mass. The force res- 
olution of the simulations is too poor to resolve convinc- 
ingly the dark matter halos that would be associated with 
individual galaxies within the cluster (one grid cell in the 
particle-mesh calculation is of order 45 kpc in length) and, 
so, it is not possible to calculate the line of sight velocity 
dispersion of member galaxies directly. However, in the 
absence of significant velocity bias in both observed and 
high-resolution numerical clusters (e.g. Lubin & Bahcall 
1993; Bromley ct al. 1995; Ghigna et al. 1998), a random 
subset of the particles can be drawn from each cluster to 
estimate the velocity dispersion that would be expected 
for the member galaxies. 

Each cluster was viewed from 1000 random orientations 
and the line of sight velocity dispersion, a v , was computed 
as a function of projected radius relative to the cluster 
center of mass. Two types of annuli were used for the 
computation: independent annuli (i.e. a v was computed 
in thin annuli with differential radius, R) and cumulative 
annuli (i.e. o~ v was computed in wide annuli which shared 
a fixed inner radius, i? m in, and differed only by the maxi- 
mum radius of the annuli, i? m ax)- That is, the use of the 
independent annuli yields a measurement of a v at a par- 
ticular projected distance from the cluster center while the 
use of the cumulative annuli yields a measurement of o~ v 
averaged over the entire cluster (out to some maximum ra- 
dius). Throughout, the minimum radius from the cluster 
centers of mass, i? m in, was taken to be a distance equal to 
the length of two grid cells in the particle-mesh simulation 
since below that scale the gravitational force is softened 
by the N-body computational technique. 

Shown in Figs. 6 and 7 (crosses) are the mean values of 
g v that were calculated directly from the line of sight ve- 
locities of particles within the clusters. Independent annuli 
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were used in Fig. 6 and cumulative annuli were used in Fig. 
7. The error bars in the figure show the formal l-a disper- 
sion amongst the different projections of the clusters. The 
velocity dispersion computed using independent annuli de- 
creases monotonically with radius in all three clusters but 
the decrease is slow enough such that averaged over large 
scales within the clusters (i.e. a v computed in the cumu- 
lative annuli) the velocity dispersion is roughly constant. 
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Fig. 6: Line of sight cluster velocity dispersions, a v , as function of pro- 
jected radius. Crosses show the mean value of <7 V computed directly from 
the velocities of random subsets of the constituent particles and the er- 
ror bars show the formal \-a deviation amongst 1000 random lines of 
sight. Squares show the value of a v inferred for the clusters on the basis 
of the mean weak lensing shear, under the assumption that the cluster 
potential is well-represented by an isothermal sphere; error bars show 
the formal \-a deviation amongst the 10 random lines of sight for which 
direct ray tracing was performed. In this figure a v has been computed 
using independent annuli with proper radius R. 

Also shown in Figs. 6 and 7 (squares) are the mean val- 
ues of a v that are inferred for the clusters on the basis of 
the weak lensing shear field, assuming that the cluster po- 
tentials can be well-represented by isothermal spheres (e.g. 
Fig. 5). From the 10 different projections for which direct 
ray tracing was performed, the mean shear was computed 
using both independent and cumulative annuli identical 
to the annuli used to compute the velocity dispersions of 
the particles themselves. The values of "f(R) and 7(i? ma x) 
obtained from the ray trace analysis were then used in con- 
junction with equations (9) and (10) to infer the variation 
of the cluster velocity dispersion with radius. Error bars 
in Figs. 6 and 7 show the formal l-a dispersion amongst 
the different cluster projections. In contrast to the ve- 
locity dispersion measured directly for the particles, the 
velocity dispersion inferred from the weak lensing analysis 
increases with radius monotonically. 
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Fig. 7: Same as Fig. 6 except that in this figure a v has been computed 
using large cumulative radii of outer radius i? max (see text). 



3.3. Cluster Mass Estimates 

Here we compute mass profiles for the clusters using 
the following simple estimators: (1) the mean value of the 
weak lensing shear under the assumption of an isothermal 
cluster potential, (2) the dynamical mass obtained from 
the line of sight velocity dispersion of the particles under 
the assumption of an isothermal cluster potential, and (3) 
the classical virial estimator. The cluster mass profiles ob- 
tained using the estimators are compared directly to the 
true mass profiles (e.g. Fig. 4) and throughout we will plot 
ratios of the estimated and true cluster mass profiles as a 
function of radius. 

Shown in Figs. 8 and 9 are the mass profiles obtained 
from the mean weak lensing shear under the assump- 
tion of an isothermal cluster potential. Results for the 
2-dimensional projected mass profiles arc shown in Fig. 
8 and the 3-dimensional mass profiles are shown in Fig. 
9. The velocity dispersion, a v (R), inferred from the 
circularly-averaged weak lensing signal (e.g. the squares 
in Fig. 6 and 7) was used in equations (7) and (8) above 
to compute M(r)i ens and M(R)\ CDS . In both Figs. 8 and 9 
the circles indicate that the value of a v (R) was determined 
using the large, cumulative annuli (i.e. an average velocity 
dispersion over the cluster out to a maximum radius of R) . 
The solid squares in these figures indicate that the value 
of <J V (R) was determined using thin, independent annuli 
(i.e. a value of the velocity dispersion computed at a par- 
ticular distance, R, from the cluster center of mass). From 
the weak lensing analysis it is not possible to measure the 
direct dependence of a v on the 3-dimensional radius, r, 
and, so, to compute the 3-dimensional mass profile we have 
taken the velocity dispersion to be a v (r) = <r v (R = r). 
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Fig. 8: The 2-dimcnsional, projected cluster mass profile obtained from s~ 

the weak lensing analysis compared to the true cluster mass profile. Solid 3 

s 

squares indicate that the value of a v used in equation (8) was determined ^ 
from independent annuli of differential radius R\ open circles indicate 
that the value of a v was determined from large, cumulative annuli with 
outer radii of i? m ax — R- The data points shown by the squares have 
been plotted such that R is the value of the projected radius at the mid- 
points of the independent radial bins and the data points shown by the 
circles arc plotted such that R is the value of i? max (i-c. for the circles 
R corresponds to the outermost radius of the annulus used in the cal- 
culation). Error bars show the l-<7 dispersion in M(R) amongst the 10 
different projections for which ray tracing was performed. 

Fig. 8 shows that there is clearly a scale-dependent sys- 
tematic deviation of the 2-dimcnsional projected cluster 
mass profile determined from the simple weak lensing anal- 
ysis adopted here. Overall the trend is for M(R)\ CDS to 
increase monotonically with radius, underestimating the 
true projected mass at small radii and overestimating the 
true projected mass at large radii. The overestimate of the 
projected mass at large radii is simply a reflection of the 
fact that the isothermal sphere is, by definition, infinite 
in extent while the actual clusters are confined to a finite 
radius of f2oo- (Note, however, that we performed a few 
test cases in which the proper radius of the numerical clus- 
ters was increased to a value of r = 4 Mpc and this had a 
negligible effect upon the measured shear and, hence, the 
inferred projected mass.) 

In contrast to the results for M(R)\ cns , there is only a 
weak scale dependence in the deviation of M (r)i ons from 
the true 3-dimensional mass. Over most scales there is 
quite good agreement between the true cluster mass pro- 
files and M(r)i ons as determined from values of cr v that 
were computed using independent annuli. When values of 
<j v determined from the large cumulative annuli are used, 
M{r)\ ens systematically underestimates the true cluster 



mass on scales significantly less than r2oo- At large radii, 
however, M(r)\ ens is in very good agreement with the true 
mass of the cluster for the case in which a v was determined 
from the large, cumulative annuli (i.e. a v is determined 
from the mean shear over the entire cluster). This result 
may seem a bit surprising given the fact that the clusters 
are better represented by NFW density profiles than they 
are by isothermal spheres. However, for NFW-type ob- 
jects with masses comparable to those of our numerical 
clusters, the mean shear interior to -R200 differs relatively 
little (< 10%) from that of an isothermal sphere that has 
an identical mass contained within r2oo (Oaxaca Wright 
& Brainerd 1999). Hence, the isothermal sphere approx- 
imation should yield a reasonable estimate of the cluster 
mass contained within r2oo, provided the mean shear is 
computed interior to a projected radius of i?200- 
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Fig. 9: The 3-dimcnsional cluster mass profile obtained from the weak 
lensing analysis compared to the true cluster mass profile. Solid squares 
indicate that the value of <r v used in equation (7) was determined from 
independent annuli of differential radius R; open circles indicate that the 
value of (7 v was determined from large, cumulative annuli with i? max — T. 
The data points shown by the squares have been plotted such that r is 
the value of the 3-dimensional radius at the midpoints of the indepen- 
dent radial bins and the data points shown by the circles are plotted 
such that r is the value of r — i? max . Error bars show the 1-er dispersion 
in M(r) amongst the 10 different projections for which ray tracing was 
performed. 

Shown in Figs. 10 and 11 are the cluster mass profiles 
obtained from equations (7) and (8) in which a v is taken to 
be the mean particle velocity dispersion measured directly 
from random subsets of particles. The 2-dimensional pro- 
jected mass profile, M av (R), is shown in Fig. 10 and the 
3-dimensional mass profile, M av (r), is shown in Fig. 11. 
As in Figs. 8 and 9, circles refer to values of a v computed 
using the large cumulative annuli and squares refer to val- 
ues of <r v computed using the thin, independent annuli. 
Both the projected mass profiles and the 3-dimensional 
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mass profiles estimated directly from the particle velocity 
dispersions show scale-dependent deviations from the true 
mass profile. In this case the cluster mass is overestimated 
at very small radii and underestimated over most of the 
cluster. 



\ t * ( ) 
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Fig. 10: The 2-dimcnsional, projected cluster mass profile obtained di- 
rectly from the measured particle velocity dispersion (assuming the clus- 
ters to be isothermal spheres) compared to the true cluster mass profile. 
Solid squares indicate that the value of a v used in equation (8) was de- 
termined from independent annuli of differential radius R\ open circles 
indicate that the value of a v was determined from large, cumulative an- 
nuli with outer radii of i? max — R. Error bars show the l-cr dispersion in 
M(R) amongst the 1000 projections from which the mean linc-of-sight 
velocity dispersion was computed. 

Lastly, shown in Fig. 12 is a 3-dimensional mass profile 
computed for the clusters using a virial mass estimator. 
The classical cluster virial mass estimator is: 



M = 



2G 



(13) 



where R e is the mean effective radius as projected on the 
sky: 

1 N 1 



\Ri - Rj\ 



and N is the number of galaxies in the cluster. Again, we 
cannot resolve the individual dark matter halos of mem- 
ber galaxies and, so, the virial analysis was performed on 
the clusters using random subsets of the particles. Parti- 
cles contained within concentric spheres of radius r cen- 
tered on the cluster centers of mass were viewed from 1000 
random orientations and a v and R^ 1 were computed for 
each orientation. Values of r were increased incremen- 
tally to r max = r 2 oo, and M(r), the mass contained within 
the concentric spheres, was computed using equation (13) 



above. From Fig. 12, the virial mass estimator leads to 
a scale-dependent deviation from the true 3-dimensional 
mass profile in the sense that the cluster mass is overesti- 
mated on small scales. On large scales (and, in particular, 
near the "edges" of the clusters), however, the virial mass 
estimator reproduces the true cluster mass quite well. 



r [Mpc] 

Fig. 11: The 3-dimensional cluster mass profile obtained directly from 
the measured particle velocity dispersion (assuming the clusters to be 
isothermal spheres) compared to the true cluster mass profile. Solid 
squares indicate that the value of a v used in equation (7) was deter- 
mined from independent annuli of differential radius R; open circles in- 
dicate that the value of a v was determined from large, cumulative annuli 
with outer radii of i? max — R. Error bars show the 1-a dispersion in 
M(r) amongst the 1000 projections from which the mean linc-of-sight 
velocity dispersion was computed. 

It should be noted that gravitational force softening in 
the numerical simulation will, necessarily, affect dynami- 
cal mass estimates of simulated objects (see, e.g., Tormen, 
Bouchet & White 1997). That is, on scales smaller than 
or of order the smoothing length, the mass will be severely 
overestimated simply due to numerical effects. We have, 
therefore, restricted our analyses to radii at which the ef- 
fects of force softening on the mass estimate should be 
small. In particular, any overestimate of the mass caused 
by numerical effects is expected to be at most of order 3% 
to 4% in the innermost radial bins and will drop rapidly 
to zero for the bins with larger radii. 

Cen (1997) and Reblinsky & Bartelmann (1999) have 
also investigated the virial masses obtained for numerical 
clusters, though not for objects as massive as those pre- 
sented here. Reblinsky & Bartelmann (1999) find that the 
virial mass severely overestimates the true masses of clus- 
ters whose masses are less than a few times 10 14 M Q . The 
degree of overestimation decreases with increasing clus- 
ter mass, however, and appears to converge in the mean 
to the true cluster mass for their most massive objects. 
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Reblinsky & Bartelmann's results are broadly consistent 
with those of Cen (1997), though differences in the proce- 
dures used to select and analyse the clusters makes direct 
comparisons between the two not entirely straightfoward. 
Direct comparisons between our results and those of Cen 
(1997) and Reblinsky & Bartclmann (1999) are also not 
straightforward. 
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Fig. 12: The 3-dimcnsional cluster mass profile obtained from the clas- 
sical virial mass estimator compared to the true cluster mass profile. 
Particles contained within concentric spheres of radius r centered on the 
cluster center of mass were used to determine the mean values of R e and 
(T„ required for the evaluation of equation (13). Error bars show the 1-c 
dispersion in M(r) amongst 1000 random projections. 

In our analyses above we have expressly calculated R e 
for each of the subsets of the particles, whereas Cen (1997) 
and Rcblinksy & Bartclmann (1999) do not. Also, we 
have selected our clusters from a 3-dimensional mass dis- 
tribution while Cen (1997) and Reblinsky & Bartelmann 
(1999) select their clusters based on 2-dimensional projec- 
tions and the assignment of luminous galaxies to a random 
subset of the particles in their simulations. As such, their 
analyses attempt to at least partially address the issue of 
contamination by interloper galaxies and false detections 
of clusters in the limit of realistic observational data. In 
contrast, our results above are effectively derived in the 
limit of ideal data (i.e. the values of R e and al are com- 
puted from objects which are known a priori to be con- 
tained within the cluster under investigation). 

4. DISCUSSION 

The cluster mass results which are the most relevant 
for direct comparison to observational investigations are 
those that were obtained using large, cumulative annuli 
(i.e. the shear and particle velocity dispersion averaged 
over large scales in the cluster) as well as the virial es- 
timate in which R e is the effective radius determined for 



the "entire" cluster. Although in principle the shear and 
velocity dispersion can be measured at independent radii 
in observed clusters, the data are generally too sparse and 
noisy for this to be practicable. (See, however, Bonnet et 
al. (1994), Tyson & Fischer (1995), Squires et al. (1996b), 
Fisher & Tyson (1997) and Carlberg et al. (1997b) for ex- 
ceptions to this.) 

The mass profiles plotted in Figs. 8 through 12 extend 
to a maximum cluster radius equal to r 2 oo (or -R200 m the 
case of the projected mass profiles). In all cases M(r) and 
M{R) in these figures refer to the mass contained within a 
3-dimensional radius, r, or a projected radius, R. Since we 
have defined the clusters to consist of all particles inside 
of a radius r 2 oo, we will define the total mass of a clus- 
ter to be the mass contained within this radius, M(j2oo)- 
The total mass obtained for each cluster from each of the 
estimators is, therefore, indicated by the points in Figs. 8 
through 12 that are plotted at the largest occurring values 
of the radius. 

In terms of estimating the total cluster mass (i.e. the 
mass of the cluster contained within a 3-dimensional radius 
°f r 20o), the classical virial estimator is found to be very 
successful. The total mass of the cluster is systematically 
underestimated, but only by ~ 10%. This result is some- 
what surprizing given the fact that within r2oo the cluster 
mass distributions are not perfectly smooth and substruc- 
ture exists at a significant level. Additionally, moment 
of inertia analyses performed using all particles within a 
radius r2oo of the cluster centers of mass show the clus- 
ter mass distributions to be clearly triaxial, rather than 
spherical (see Brainerd, Goldberg & Villumsen 1998 for 
the relevant discussions). Our result, therefore, suggests 
that at least in the limit of ideal data the classical virial 
mass estimator is quite robust to modest deviations from 
pure spherical symmetry and the presence of substructure 
within a cluster. 

The "isothermal" dynamical mass estimate, in which the 
measured line of sight velocity dispersion is used to infer 
the mass under the assumption of an isothermal potential, 
yields a poor estimate of the total mass of the cluster. The 
value of M(r2oo) is underestimated by ~ 40% for the case 
in which a v is determined from an average over the entire 
cluster and is underestimated by ~ 70% for the case in 
which a v is computed at a projected radius of R = i?2oo 
(i.e. Fig. 11). 

Provided the mean shear used to infer the cluster veloc- 
ity dispersion is computed using a large, cumulative annu- 
lus in which the shear is averaged over the entire cluster, 
the weak lensing estimate of M(r 2 oo) is found to be in 
excellent agreement with the total cluster mass (i.e. Fig. 
9, open circles). In contrast, however, the shear measured 
solely at a radius of R — i?2oo yields a ~ 25% overestimate 
of the total cluster mass (i.e. Fig. 9, solid squares) due to 
the fact that the clusters are finite in extent, rather than 
infinite. 

Because of its promise to yield direct measurements 
of the masses of galaxy clusters independent of dynam- 
ics and hydrodynamics, weak lensing mass estimates of 
cluster masses are currently of particular interest. Given 
the fact that most high-quality observations of the weak 
lensing shear due to clusters have been obtained only on 
relatively small scales (i.e. radii significantly less than 1 
or 2 Mpc), Figs. 8 and 9 suggest some caution regarding 
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the interpretation of recent cluster mass estimates that are 
based on a measurement of an average value of the weak 
shear together with an assumption of an isothermal clus- 
ter potential. In particular, a measurement of the mean 
shear in which the mean is computed within an aperture 
whose outer radius is significantly less than -R200 yields 
a mass estimate that differs systematically from the true 
mass. For example, the projected mass within a radius of 
0.5 Mpc, M(R = 0.5 Mpc), is underestimated by ~ 40% 
and the 3-dimensional contained mass, M(r = 0.5 Mpc), 
is underestimated by ~ 35%. Interestingly, in an analy- 
sis of observed cluster lensing data, Wu et al. (1998) found 
that weak lensing mass estimates that were performed over 
small cluster radii did seem to underestimate the contained 
mass in a systematic manner. 

The results for M\ ens (r) shown in Fig. 9 are, however, 
encouraging at large cluster radii. In particular, with the 
advent of large format CCD cameras capable of wide-field 
imaging, it will be possible to measure the weak shear due 
to lensing clusters at radii of order a few Mpc in a reason- 
ably routine fashion. An example of such deep wide-field 
imaging is the data obtained with the UH 8K CCD mosaic 
camera which has recently resulted in a detection of large 
scale coherent weak shear in the images of ~ 30, 000 faint 
background galaxies due to lensing by the supercluster 
MS0302+17 (Kaiser et al. 1998). Given the apparent uni- 



versality of the Navarro, Frenk & White density profile (i.e. 
dissipationless collapse generically leads to the formation 
of an object with an NFW profile), our results suggest that 
it will be possible to estimate a total 3-dimensional clus- 
ter mass fairly accurately with wide-field imaging simply 
by computing the mean of the shear over the entire clus- 
ter and adopting an isothermal lens potential. In the short 
term, such observations will hopefully provide a resolution 
to the remaining discrepancies between cluster masses es- 
timated from weak lensing and virial techniques. (This is, 
of course, providing that the redshift distribution of the 
lensed galaxies is well-constrained and is not, in itself, a 
large source of uncertainty in the interpretation of the ob- 
served shear.) In the long term, large surveys from which 
the weak lensing shear can be detected out to large cluster 
radii should have the ability to yield uniform samples of 
objects, including a reasonably accurate mass-selection cri- 
terion, without necessarily requiring a full reconstruction 
of the density profile of each individual lensing cluster. 
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